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The derivative discontinuity of the exchange-correlation functional of density-functional theory 
is cast as the difference of two types of electron affinities. We show that standard Kohn-Sham 
calculations can be used to calculate both affinities, and that their difference benefits from substantial 
and systematic error cancellations, permitting reliable estimates of the derivative discontinuity. 
Numerical calculations for atoms and molecules show that the discontinuity is quite large (typically 
several eV) , and significantly and consistently improves the agreement of the calculated fundamental 
energy gap with the experimental gaps. The percentage error of the single-particle gap is reduced 
in favorable cases by more than one order of magnitude, and even in unfavorable cases by about a 
factor of two. 
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From spectroscopy to transport, there is hardly any 
property of a quantum many-particle system that does 
not in some way depend on whether there is a gap in 
the energy spectrum, and what the size of this gap is. 
The fundamental gap is a ground-state property of the 
TV-body system, defined [j], 0, 0] in terms of the ground- 
state energy E(N) as E g = [E(N - 1) - E(N)} - [E(N) - 
E(N + 1)], where E(N - 1) - E(N) is the energy change 
upon removing the iV'th particle from the N particle 
system and E(N) — E(N + 1) that upon removing the 
N + 1st particle from the N + 1 particle system. In 
a noninteracting system, this definition reduces to the 
familiar energy gap between single-particle levels. 

For interacting systems, approximate many-body cal- 
culations of total energies and energy gaps are typ- 
ically performed within the framework of density- 
functional theory (DFT) [l|, 0, [I| or, for small finite sys- 
tems, Hartree-Fock (HF) and post-Hartree-Fock meth- 
ods. DFT provides, in addition to the ground-state den- 
sity and related quantities, also a set of single-particle 
eigenvalues, the so-called Kohn-Sham (KS) spectrum. 
The difference between the energy of the highest occu- 
pied and lowest unoccupied of these single-particle levels 
is the KS gap, which in extended systems become the 
band structure gap. Neither the KS nor the HF single- 
particle gaps correspond to the experimental gap, the 
former typically being too small and the latter too large. 

Generally, one can write E g = E^ + A xc , which de- 
fines A xc as the difference between the exact fundamen- 
tal gap and the exact KS single-particle gap. In atomic 
physics and quantum chemistry, the importance of a 
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nonzero A xc for chemical hardness is well known 0, IH [f| . 
Neglect of A xc has also been shown to lead to large errors 
in the calculation of Rydberg excitations In H , charge- 
transfer processes and transport [1, ITol. and the 
ionization probability of atoms exposed to strong electro- 
magnetic fields 12|. In semiconductors, approximate en- 
ergy gaps calculated in DFT often drastically underesti- 
mate the experimental gap [HI, EH- In Mott insulators, in 
particular, the entire gap is due to A^c [l6[. Recently 
it was pointed out that a similar discontinuity problem 
can also appear in the spin-DFT calculation of spin gaps 
in materials relevant for spintronics [T?| . The question 
whether the neglect of A xc or the error in j^KS,approx 
is responsible for the underestimate of the band gap in 
solids is considered in a standard textbook in the field 
to be 'of tremendous practical importance' and the 
calculation of A xc is ranked in a recent monograph as 
' certainly one of the major outstanding problems in con- 
temporary DFT [3[, but no general answer is known. 

In the present paper we draw attention to an alterna- 
tive representation of A XCl which casts it as a difference 
of single-particle eigenvalues, similar to the KS gap. We 
point out that this relation provides a simple physical in- 
terpretation of the elusive xc discontinuity, and use it to 
estimate xc discontinuities of atoms and molecules. The 
resulting correction A^ s * substantially improves agree- 
ment with experimental fundamental gaps, reducing the 
percentage deviation from experiment by more than a 
factor of 10 in favorable cases and by about a factor of 
two in less favorable ones. 

The standard representation of A xc is based on ensem- 
ble DFT for open systems, where all three quantities in 
Eg = Eg <s + A xc can be related to derivative disconti- 
nuities of universal density functional [H, El, Eil ■ The 
fundamental gap is the derivative discontinuity of the to- 
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TABLE I: Experimental gap, Eg Xp = I exp — A exp , KS single-particle gap, Ef , and KS gap corrected by adding our estimate 
for the discontinuity, Eg FT = Eg + A^*, for atoms Li (Z = 3) to Ca (Z = 20), with exception of the nobel gases. The 
values for A% s ^ were obtained from Eq. ([7Jl using the B88LYP functional and the 6-311G(d) basis sets. All data in eV. 
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In these equations rj stands for an infinitesimal variation 
of the system's particle number. Equation (|3|) shows that 
A xc is a system-dependent shift of the xc potential v xc (r) 
as it passes from the electron-poor to the electron-rich 
side of integer N. 

Theses three relations are useless to calculate gaps 
from most currently available approximate density func- 
tionals, which typically have no discontinuities. For two 
of the three quantities above, alternative ways of calcu- 
lation, more useful in practice, are widely known. Total 
energies are calculated easily from DFT, so that one can 
employ the definition of ionization energy I = E(N — 
1) - E(N) and electron affinity A = E(N) - E(N + 1) 
to calculate the fundamental gap from 

E g = E(N + 1) +E(N- 1) - 2E(N) = I- A. (4) 

Single-particle energies are obtained from the KS equa- 
tion, as a byproduct of calculating the total energies, and 
yield the KS single-particle gap 



E« s = e N+1 (N) 



(5) 



where cm(N) denotes the M'th eigenvalue of the N- 
electron system. The third term, the xc discontinuity 
A xc , has resisted all attempts of describing it directly 
by common density functional, such as LDA and GGA, 
which are continuous as a function of N and thus have 
no xc discontinuity. 

However, we note that / and A can be calculated in 
DFT not only from ground-state energies, but also from 



single-particle eigenvalues, via / = — ejv(iV) and A — 
—cn+i(N + 1) (the analogue of Koopmans' theorem in 
DFT) G], [2(3] . By using these relations and the definition 
of E^ one finds, upon combining Eqs. (j4]) and ([5]) [HI], 



E„ 



e N+1 {N)-e N (N) + A x 



e N+1 (N +l)-e N (N), 
(6) 

which implies 

A xc = e N+1 (N + l)-e N+1 (N) = A KS -A. (7) 

In the last step we used the fact that the affinity of the 
KS system, Aks> is simply the negative of the energy of 
the lowest unoccupied orbital. 

One way to estimate A xc is by subtracting Eq. ([5]) from 
(|4]). In the following, we use Eq. (JTJjwhich is equivalent 
in principle, but simpler in practice [221 ]. It also provides 
an intuitive interpretation of the discontinuity: in an in- 
teracting system, the electrons repell, hence the energy 
cost of removing the outermost electron from the nega- 
tive species (which is measured by the electron affinity) 
is reduced, A < AkSi an d a positive A xc results. 

If the right-hand side of Eq. (O could be calculated 
exactly, this procedure would determine the exact xc dis- 
continuity. An estimate of A xc is thus obtained by using 
in ||7J| the KS eigenvalues obtained in two approximate 
KS calculations, one for the neutral species, the other for 
the anion. 

Such approximate calculations involve two distinct 
types of errors, one associated with the approximations 
used for the xc functional, the other with the finite size 
of the basis set. As a consequence, each of the two affini- 
ties in Eq. ([7]) is predicted wrongly by standard com- 
binations of functionals and basis sets. Typically, the 
self-interaction error inherent in common LDA and GGA 
type functionals shifts the eigenvalues up, in some cases 
so much that the anions become unbound [23[. On the 
other hand, the finiteness of the basis set artificially sta- 
bilizes the anion (23|. As a consequence of this error 
cancellation, practical methods for calculating affinities 
from LDA and GGA are available [H, [H [H [2^. 

Our key argument is that derivative discontinuities are 
protected from functional errors and basis-set errors by 
a distinct additional error compensation, independent of 
the one just described. Namely, Eq. (0) casts A xc as 
a difference of two affinities. Even if each is predicted 
wrongly on its own, their difference is expected to bene- 
fit from substantial additional error cancellation. In fact, 
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if all KS eigenvalues are shifted up by roughly the same 
amount, energy differences are preserved, and even posi- 
tive eigenvalues (unbound anions) can provide reasonable 
discontinuities for the bound neutral system. 

We call this calculation of A xc by means of Eq. ([7|) 
an estimate, and not an approximation, to stress that it 
exploits an error cancellation that is hard to quantify a 
priori. However, in recent work on models of harmon- 
ically confined systems [27j this estimate was found to 
lead to significantly improved gaps. Here we explore the 
performance of Eq. ([Jj in ab initio calculations for atoms 
and molecules. 

In Table H] we compare, for 16 light atoms, the experi- 
mental gap, the KS single-particle gap, and the DFT gap, 
defined as the sum of the KS gap and the estimated xc 
discontinuity. The KS calculations were performed with 
the GAUSSIAN 98 |H| program, using the B88-LYP func- 
tional and the 6 — 311G(d, p) basis sets. Table Q] shows 
that the error of the KS gap is significantly and consis- 
tently reduced by adding the estimated xc discontinuity 
to the KS gap, dropping by more than an order of mag- 
nitude — from 7.67eV to 0.606el/ — over the data set 
in Table |TJ This large drop, together with the fact that 
the improvement is systematic (obtained not only on av- 
erage, but in every individual case), strongly suggests 
that Eq. ([7]) is a reliable and useful way of obtaining the 
discontinuity. 

Figure [1] is a plot of the data in Table [I] revealing 
that A xc roughly follows the atomic shell structure. Par- 
ticularly small discontinuities are found for atoms with 
one electron outside a closed shell, such as Li and Na. 
The largest discontinuities are, however, not observed for 
closed-shell systems but for systems one electron short of 
a closed shell, as is seen comparing F with iVe or CI 
with Ar. We interpret this by means of Eq. as a 
consequence of the fact that A xc is related to two affini- 
ties, which involve negative species with one additional 
electron, leading to a closed shell for F ~~ and Cl~ . 

Specifically for the Be atom, we can further compare 
with independent theoretical expectations, because the 
discontinuity of this atom has previously been estimated 
by Jones and Gunnarsson (JG) [29] by comparing the 
experimental gap to a near-exact KS gap obtained ear- 
lier by Pedroza and Almbladh [3(| from CI densities and 
approximate inversion of the KS equation. Our value 
Af c e = b.heV is encouragingly close to Af c e ' JG = 5.7eV. 

Next, we turn to molecules. In Table Ql] we compare 
our estimate of A xc to many-body values of A xc and 
to experimental fundamental gaps. The many-body dis- 
continuity is obtained (3~0 ] by performing coupled-cluster 
calculations to generate a near-exact density, followed by 
inversion of the KS equation to obtain the correspond- 
ing near-exact KS potential, solution of the KS equation 
with that potential to obtain the near-exact KS gap, and 
subtraction of that gap from the experimental fundamen- 
tal gap [3l[. The first step is impractical for larger sys- 
tems, whereas the last step involves using the experimen- 
tal gap, which makes the method empirical. For these 
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FIG. 1: Plot of the data in TablelU Kohn- Sham single-particle 
gap (triangles), experimental fundamental gap (open squares) 
and Kohn-Sham single-particle gap corrected by adding our 
estimate of the discontinuity (full circles), for 16 light atoms. 
The lines are guides for the eye, illustrating that A% s c re- 
covers not only the overall value, but also fine details of the 
behaviour as a function of Z that were lost in the KS gap. 



reasons, the simple estimate obtained from Eq. (|7|) may 
constitute a useful alternative, provided it turns out to 
be sufficiently reliable. 

In fact, the estimated value of A e x s * depends on the cho- 
sen basis set and functional, as well as on whether the 
anion geometry is separately optimized (leading to adia- 
batic affinities and discontinuities) or hold fixed at that 
of the neutral species (vertical affinities and discontinu- 
ities). Tests of different combinations of methodologies 
indicate that best (i.e., closest to experiment) gaps are 
obtained if the discontinuity is calculated from vertical 
affinities. 

Table HI] shows that the systems fall in two classes. For 
one class, comprising CO, H 2 CO, H 2 S, HCN, N 2 , PH 3 , 
Cl 2 , S0 2 , and C2H4, the gap error with respect to exper- 
iment is reduced by a similar margin as for atoms, or even 
more, and estimated and calculated A xc agree well. For 
the other class, comprising H2O, NH%, HF, CH4 and 
perhaps C2H2, the percentage error of the gap drops by 
a factor of two, instead of by one order of magnitude, and 
the estimate recovers about 50% of the many-body value 
of A xc . An empirical indicator of which class a system 
belongs to is the sign of the KS affinity: if this is negative 
(i.e., the lowest unoccupied KS orbital has positive eigen- 
value) the system belongs to class II. We note that this is 
not a stability criterium because it employs the (unphys- 
ical) KS affinity and because for all molecules in Table HT1 
(with exception of CI2 and S0 2 ) the anionic species is 
experimentally unstable. Rather, it indicates a partial 
loss of the error cancellation on which our use of Eq. ([7]) 
is based. Nevertheless, even in these "unfavorable" cases, 
the estimate still provides a systematic improvement on 
the KS gap, reducing the error with respect to experi- 
ment by a factor of two. For both classes, A x s * clearly 
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TABLE II: Comparison of calculated and experimental gaps for small molecules. First column: experimental gap. Second 
column: KS single-particle gap. Third column: percentage deviation of KS gap from experimental gap. Fourth column: 
derivative discontinuity estimated from single-particle eigenvalues, obtained from the B88-LYP functional on the 6-311G(d,p) 
basis sets. Fifth column: DFT gap. Sixth column: percentage deviation of DFT gap from experimental gap. Seventh column: 
many-body estimate of the derivative discontinuity. All values in eV. 
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provides a quantitative correction to single-particle gaps, 
which may be useful in improving, e.g., the DFT descrip- 
tion of chemical hardness 0, Q or of transport through 
single molecules [1, 0] . 

In summary, we have cast the derivative discontinuity 
of DFT as a difference of two affinities, Eq. (JT]), and used 
approximate KS calculations of these to obtain estimates 



for the discontinuity in atoms and molecules. Our results 
are consistent with previous results, where available, and 
significantly and consistently reduce the error between 
calculated and measured fundamental gaps. 

This work was supported by FAPESP, CAPES and 
CNPq. 



[1] R. M. Dreizler and E. K. U. Gross, Density Functional 

Theory ( Springer- Ver lag, Berlin, 1990). 
[2] R. G. Parr and W. Yang, Density-Functional Theory of 

Atoms and Molecules (Oxford University Press, Oxford, 

1989). 

[3] G. F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Liquid (Cambridge University Press, 2005). 

[4] A. Gorling and M. Levy, Phys. Rev. A 52, 4493 (1995). 

[5] D. J. Tozer and F. De Proft, J. Phys. Chem. A 109, 8923 
(2005). 

[6] A. Wasserman, N. T. Maitra and K. Burke, Phys. Rev. 

Lett. 91, 263001 (2003). 
[7] D. J. Tozer and N. C. Handy, J. Chem. Phys. 109, 10180 

(1998). 

[8] C. Toher, A. Filippetti, S. Sanvito and K. Burke, Phys. 

Rev. Lett. 95, 146402 (2005). 
[9] S.-H. Ke, H. U. Baranger and W. Yang, J. Chem. Phys. 
126, 201102 (2007). 
[10] D. J. Tozer, J. Chem. Phys. 119, 12697 (2003). 
[11] Q. Wu and T.VanVoorhis, Phys. Rev. A 72, 024502 
(2005). 

[12] M. Lein and S. Kiimmel, Phys. Rev. Lett. 94, 143003 
(2005). 

[13] R. W. Godby, M. Schliiter and L. J. Sham, Phys. Rev. 
Lett. 56, 2415 (1986). 



[14] M. Griming, A. Marini and A. Rubio, J. Chem. Phys. 

124, 154108 (2006). 
[15] J. P. Perdew and M. Levy, Phys. Rev. Lett. 51, 1884 

(1983). 

[16] N. A. Lima, L. N. Oliveira, and K. Capelle, Europhys. 

Lett. 60, 601 (2002). 
[17] K. Capelle and G. Vignale, Phys. Rev. Lett. 86, 5546 

(2001). 

[18] L. J. Sham and M. Schliiter, Phys. Rev. Lett. 51, 1888 
(1983). 

[19] W. Kohn, Phys. Rev. B 33, 4331 (1986). 
[20] Q. Wu, P. W. Ayers, and W. Yang, J. Chem. Phys. 119, 
2978 (2003). 

[21] L. J. Sham and M. Schliiter, Phys. Rev. B 32, 3883 
(1985). 

[22] Eq. (J7J) does not require calculations for the iV — 1 particle 
system, and also avoids numerically troublesome calcu- 
lation of A xc as difference of much larger numbers. 

[23] N. Rosch and S. B. Trickey, J. Chem. Phys. 106, 8940 
(1997). 

[24] J. M. Galbraith and H. F. Schaefer, III, J. Chem. Phys. 

105, 862 (1996). 
[25] G. de Oliveira, J. M. L. Martin, F. de Proft, and P. Geer- 

lings, Phys. Rev. A 60, 1034 (1999). 
[26] O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 122, 



■5 



184107 (2005). 

[27] K. Capelle, M. Borgh, K. Karkkainen and S. M. 
Reimann, Phys. Rev. Lett. 99, 010402 (2007). 

[28] Gaussian 98, M. J. Frisch et al., Gaussian Inc., Pitt- 
burgh, PA (1998). 

[29] R. O. Jones and O. Gunnarsson, Rev. Mod. Phys. 61, 



689 (1989). 

[30] C.-O. Almbladh and A. C. Pedroza, Phys. Rev. A 29, 
2322 (1984). 

[31] M. J. Allen and D. J. Tozer, Mol. Phys. 100, 433 (2002). 
[32] R. G. Pearson, Inorg. Chem. 27, 734 (1988). 



